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ASSESSMENT OF A TRANSITIONAL BOUNDARY LAYER THEORY 
AT LOW HYFERSONIC MACH NUMBERS 

S. J. Shamroth and H. McDonald 
United Aircraft Research Laboratories 


SUMMARY 


This investigation assesses the accuracy of a transitional 
boundary layer theory in the low hypersonic Mach number regime. The 
theory, previously shown to be applicable to subsonic and moderately 
supersonic flows, is based upon a simultaneous solution of the integral 
turbulence kinetic energy equation and a finite -difference solution of 
the boundary layer partial differential equations of motion. The 
investigation reported here shows the capability of the original form 
of the theory to accurately predict Stanton numbers and mean velocity 
profiles through the transitional regime at low hypersonic Mach numbers 
and with a slight modification to correctly predict the effect of Mach 
number and wall cooling on transition Reynolds number. The procedure 
also shows promise of predicting initiation of transition for a given 
free-stream turbulence disturbance level; an accurate assessment of 
this capability, however, cannot be made until more detailed experimental 
data become available. An investigation of the direct effect of fluctu- 
ating pressure shows the pressure -dilitation contribution to the turbu- 
lence kinetic energy balance, usually ignored at low Mach numbers, has, 
at most, a moderate effect in the low hypersonic Mach number regime. 

In addition, transition calculations performed with the usual free- 
stream turbulence intensity source term set to zero and a magnitude 
ascribed to the turbulence source term representing the direct absorp- 
tion of 'acoustic energy by the boundary layer show that to trigger 
transition requires only a small amount ( ~ 1 %) of acoustic energy 
absorption. 


INTRODUCTION 


Transitional boundary layers play an important role in the 
successful design and operation of a hypersonic reentry vehicle. 

During transition both the wall shear and wall heating can reach their 
peak values, thus having a potentially important effect on vehicle drag 
and the amount of wall cooling required to maintain structural integrity. 
In addition, the wake structure behind the vehicle which depends upon 
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the wall boundary layer development may vary significantly with transi- 
tion location. Thus an analytical procedure capable of predicting 
transitional hypersonic boundary layers would be a useful tool in 
vehicle design. 

Presently, there exist three main approaches for developing a transi- 
tional boundary layer theory: the semiempirical approach, the stability 

theory approach, and the turbulence kinetic energy approach. Typical 
examples of the semiempirical approach for predicting transition loca- 
tion are discussed by Hairston (ref. l). Semiempirical procedures 
cannot be used with confidence to predict transitional properties for 
flow conditions significantly different from the flow conditions of the 
correlating data. This uncertainty is further heightened by the effects 
of tunnel influence upon transition experiments. As shown by Pate and 
Schueler (ref. 2) and Pate (ref. 3 ), correlations based upon transition 
data from one tunnel do not necessarily correspond to correlations 
based upon data from a different tunnel. Therefore, tunnel effects may 
severely limit the validity of semiempirical theories and make them 
inapplicable to free -flight situtations . 

In addition to the drawbacks inherent in data correlations, most 
semiempirical theories make the crude assumption that transition occurs 
instantaneously. A more realistic transition model developed by 
Harris (ref. 4 ) still requires specification of both the location and 
extent of the transition region but does not assume instantaneous 
transition. In ref. 4 Harris assumes the transitional shear stress as 
a linear combination of the laminar shear stress and a fully -developed 
turbulent shear stress, the linear combination factor being a universal 
function of the dimensionless streamwise location within the transitional 
region adapted from the low speed correlation of Dhawan and 
Narashima (ref. 5 )* This model has successfully predicted the develop- 
ment of transitional boundary layers in tne low hypersonic Mach number 
regime and represents a significant improvement over instantaneous 
transition models. However, since it heavily depends upon empiricism 
and was not intended to predict the effect of phenomena such as free- 
stream turbulence, Mach number, wall temperature, etc. upon transitional 
behavior, it does not -serve as a general boundary layer transition 
prediction procedure . 

The second possible approach which might lead to an analytical 
prediction of transitional boundary layers is based upon stability 
theory. In classical stability theory the flow is divided into a mean 
flow whose stability is the subject of the investigation and a super- 
imposed disturbance of specified form. The resulting equations are 
simplified through the boundary layer assumptions, linearized with 
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respect to the disturbance flow field, and then solved to determine 
conditions for amplification of the disturbance. Attempts to predict 
the initiation of transition from classical stability theory (for 
example, ref. 6) have had only qualitative success and this approach has 
yet to be applied successfully to the problem of transitional boundary 
layer development. The neglect of nonlinear phenomena and the assumption 
of purely two-dimensional disturbances makes classical stability theory, 
without major modification, an unlikely method of predicting the behavior 
of transitional boundary layers. 

The third approach for predicting the behavior of transitional 
boundary layers is based upon the solution of the turbulence kinetic 
energy equation. Whereas, stability theory in a sense extrapolates 
forward from laminar to turbulent flow, use of the turbulence kinetic 
energy equation is an effort to extrapolate backward from the turbulent 
to laminar regimes. Since the turbulence kinetic energy equation 
derives directly from the Navier-Stokes equations (e.g., ref. 7), with 
no further assumptions, it is fundamentally a more accurate mathematical 
description of the transition process than is the linearized two- 
dimensional disturbance equation (the Qrr -Sommer field equation) used in 
classical stability theory. Obviously, a formal relationship must exist 
between a transition theory based upon the turbulence kinetic .energy 
equation and classical stability analysis since the turbulence kinetic 
energy equation governs both the wave motion of stability theory and 
the developing turbulence . 

The turbulence kinetic energy equation has proven to be a useful 
tool in predicting the behavior of a wide variety of turbulent boundary 
layers (e.g., refs. 8 and 9) and has been extended to the transitional 
regime by Glushko (ref. 10), Donaldson, Sullivan, and Yates (ref. 11), 
and McDonald and Fish (ref. 12). The approaches of refs. 10 through 12 
are basically similar in philosophy: i.e., they integrate an equation 

(or, in the case of Donaldson, et al., a set of equations) governing 
the streamwise development of the turbulent shear stress; however, they 
differ in the relationships introduced to allow performance of the 
integration. At the present time, neither the work of Glushko (ref. 10) 
or Donaldson, et al. (ref. 11) has been developed into a practical pre- 
diction procedure although preliminary calculations by Beckwith and 
Bushnell (ref. 13) indicate that Glushko's method perhaps could be 
developed satisfactorily. The approach used by McDonald and Fish 
(ref. 12) has proven capable of accurately predicting the behavior of 
a wide variety of transitional boundary layers in the subsonic and 
low supersonic Mach number regimes. The present investigation assesses 
the ability of the theory of ref. 12 to predict the development of 
transitional boundary layers in the low hypersonic Mach number regime. 
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The existence of phenomena associated primarily -with hypersonic 
flow makes a straightforward application of any existing low Mach 
number transitional boundary layer theory uncertain in the hypersonic 
flow regime. For example, in subsonic flow the free -stream turbulent 
kinetic energy, a vorticity disturbance mode, plays an important role 
in determining transition location (refs. l4 and 15). At high Mach 
numbers, in addition to the vorticity mode, disturbances take the form 
of sound modes (ref. 2) and entropy modes. The effect of the entropy 
mode upon transition location is thought to be small except perhaps at 
very high Mach numbers; however, as shown by Wagner, Maddalon, and 
Weinstein (ref. 16), even at high Mach numbers the entropy disturbance 
may remain small. On the other hand, even at moderately supersonic 
Mach numbers the acoustic mode may play an important role in determining 
transition location. A successful hypersonic transitional boundary 
layer theory has to relate the vorticity and the acoustic modes so that 
the vorticity fluctuation arising from the acoustic disturbance could 
be incorporated into the analysis and, in addition, allow for the direct 
absorption of acoustic energy. A second example of a phenomenon 
peculiar to high Mach number flows is the direct effect of fluctuating 
pressure terms upon the turbulent kinetic energy balance. At low Mach 
numbers these quantities are usually negligible, however, in the 
hypersonic regime they may be important. In the present investigation 
the applicability of the McDonald -Fish transitional boundary layer 
theory to the low hypersonic Mach number regime is assessed. In making 
the assessment, comparisons between theory and experiment including 
velocity profiles, Stanton numbers, and integral thickness parameters 
are made and the ability of the procedure to predict the initiation of 
transition is determined. In addition, a preliminary assessment of 
the direct role of the fluctuating pressure is made. 

LIST OF SYMBOLS 

A n nonlinear coefficients, (see Eq. (A-5)) 

a n structural coefficients of turbulence 

B n nonlinear coefficients, (see Eq. (A-6)) 

Cp specific heat 

Cf skin friction coefficient 

jb sublayer damping factor 
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F 

G' 

K 

k 

L 


M 

M r 

Pr 

Pr T 

p 

Q 


Re x 

Rr 

Rr 

R 0 

St 

T 


u 

U T 


dimensionless stream function 
dimensionless temperature ratio 
constant coefficient, (see Eq. ( 39 )) 
thermal conductivity- 
dissipation length 
mixing length 

wake value of mixing length 
Mach number 

Mach number relative to free stream 

Prandtl number 

turbulent Prandtl number 

pressure 

heat flux 

turbulence kinetic energy 
radius 

Reynolds number based upon streamwise distance 

turbulence Reynolds number 

layer averaged turbulence Reynolds number 

Reynolds number based upon momentum thickness 

Stanton number 

static temperature 

total temperature 

streamwise velocity 

friction velocity 
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V transverse velocity 

w cross flow velocity 

x streamwise coordinate 

y transverse coordinate 

y + dimensionless transverse coordinate 

a indicator equal to one for axisymmetric flow, zero for 
two-dimensional flow 

(3 density ratio 

T intermittency factor 

7 ratio of specific heats 

8 boundaiy layer thickness 
8 g sublayer thickness 

8 th thermal boundary layer thickness 
8* displacement thickness 
8 * reference length 
€ turbulence dissipation 
rj dimensionless transverse coordinate 
6 momentum thickness 
/I viscosity 

i / kinematic viscosity 
Vj kinematic eddy viscosity 
p density 
T shear stress 
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<fi i , 4> 2 . 4* 3 integral functions (see Eqs. (20) through (22)) 
Subscripts 

e boundary layer edge condition 

CD free -stream condition 

w wall condition 
Superscripts 

— average quantity 
' fluctuating quantity 

THEORY 


The Basic Equations 

Within the framework of boundary layer theory, various authors 
(for example, Schubauer and Tchen (ref. 17)) have reduced the time- 
averaged Navier -Stokes equations to the compressible boundary layer 
equations of motion. For two-dimensional or axisymmetric flows, steady 
in the mean, the boundary layer approximations to the momentum, energy, 
and continuity equation become 


-- du d u d P dr 

P u dx + P v d y " ef* + dy 


( 1 ) 


p u Cp 


dT° . _ 
d\ 


r ai 
+ / 3vC p _ d7 


dy 


(Q+ u r) 


( 2 ) 


dp u r a 3pv r a 
dx + dy 


(3) 
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■where X and y are coordinates in the strearawise and transverse directions, 
U and V are velocity components in the X and y directions, p is density, 

P is pressure, Cp is specific heat, T 0 is total temperature, r is the 
radius of curvature for an axisymmetric body, and the exponent Q is zero 
for two-dimensional flows and unity for axisymmetric flows. The effec- 
tive shear stress, r, and the effective heat transfer, Q , are defined as 


T - 


_ du 
^ dy 


u v 


( 4 ) 


Q = k 


dT 

.dy 


- p Cp v‘ T ' 


(5) 


where p is viscosity, k is thermal conductivity, T is static temperature, 
In Eqs . (l) through (5) overbars indicate averaged quantities and primes 
indicate fluctuating quantities . The equations are valid for laminar, 
transitional, or turbulent flows; obviously, for laminar flows the 
primed quantities are zero. For turbulent and transitional boundary 
layers it is convenient to represent the contribution of the apparent 
turbulent stress, r T , to the total shear stress, r , by an effective 
turbulent viscosity, v r , and similarly the turbulent contribution to 
the total heat flux,Q , is represented conveniently by an effective 
turbulent conductivity, k T , such that 


du - — ; — 
P v i ~dy~ = ' P u v 


( 6 ) 


dr 

dy 


p Cp v' T 


(7) 


The turbulent conductivity, k y , is related to the turbulent viscosity, 
Vy, by a turbulent Prandtl number, Pr T , defined as 


Pr t " P ^ t 


( 8 ) 
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and the boundary layer momentum and energy equations, Eqs. (l) and (2), 
become 



In deriving Eq. (10) use has been made of the definition of total 
temperature 


T 


o 


T + 


2Cp 


( 11 ) 


With the flow laminar, Eqs. (3), (9)> and (10) are solved with i/j = 0 
to determine the boundary layer development. If the flow is transitional 
or turbulent, it is necessary to model Vj and Pr T . The specification of 
the turbulent viscosity, Vj , and the turbulent Prandtl number, Pfj , is 
carried out through the turbulence model described in detail in the 
subsequent section. 

After specification of the turbulent viscosity and Prandtl number, 
Eqs. (3), (9), and (10) are solved by first eliminating pv from the 
momentum and energy equations by application of the continuity equation. 
When the streamwise static pressure distribution, P(x), is specified, 
the momentum and energy equations, in conjunction with an equation of 
state and equations governing Vj and Pr T , form a closed set of nonlinear, 
parabolic, partial differential equations which can be solved upon 
specification of boundary conditions. The wall and free-stream boundary 
conditions employed in the solution are given by: 

at the wall, y = 0 

pv = ( pv)w, u = 0, T° : Tw or = o ( 12 ) 
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at the free -stream, y . CD 

P" ~ Pe u e 

T ° = T° (13) 

The subscripts 'w' and 'e' denote wall and free -stream conditions, 
respectively. The initial conditions for the problem are set by 
specifying the initial boundary layer displacement thickness, S* , and 
assuming the initial development is similar in the dimensionless 
coordinate, 17 , where 


V = 



(14) 


The numerical procedure used to solve the governing momentum and 
energy equations is a Hartree-Womersley approach in which streamwise 
derivatives are replaced by finite differences, the coordinate normal 
to the wall is nondimensionalized, and a stream function introduced. 

The resulting momentum equation is a third order nonlinear ordinary 
differential equation in the transverse coordinate and the energy 
equation is a second order nonlinear ordinary differential equation in 
the same transverse coordinate. At each streamwise station the non- 
linear coefficients of each equation are estimated from the solution 
at the previous station and the resulting linearized equations solved 
as two-point boundary value problems. The resulting solutions are 
used to obtain better estimates of the nonlinear coefficients and the 
entire process repeated until two successive solutions agree to within 
a specified tolerance. The chosen numerical procedure has three out- 
standing attributes vis-a-vis other available methods of solving the 
partial differential equations. First it has great flexibility in 
streamwise step size which follows from the Hartree-Womersley repre- 
sentation of the streamwise derivatives in conjunction with the choice 
of nondimensional coordinate system normal to the wall. Both very 
large and very small streamwise steps may be taken, making the procedure 
ideally suited for either laminar, transitional, or turbulent flow 
since each type of boundary layer usually demands a different choice 
of streamwise step size (in terms of boundary layer thicknesses). The 
second attribute follows from the use of an implicit finite -difference 
procedure to solve the boundary value problem normal to the wall. Not 
only does the use of the implicit procedure enlarge the range of 
streamwise step sizes which the program can handle but since a 
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diagonally dominant matrix is obtained, very fast matrix inversion 
procedures can be utilized. Furthermore, an implicit method more 
properly represents a parabolic differential equation and allows for a 
very flexible grid normal to the wall thereby enabling a large number 
of grid points to be taken in the sublayer region. It should be pointed 
out that in calculating transitional and turbulent boundary layers the 
present procedure does not fit an empirical profile between the wall and 
the edge of the sublayer but rather calculates the flow development in 
the sublayer regime from the full boundary layer equations of motion 
and, thus, can include the effects of wall bleed and roughness in a 
fundamental manner. The third attribute of the procedure arises from 
the treatment of the nonlinear terms in the equations of motion. In 
this regard it is recognized that although almost all finite -difference 
procedures represent a linearization of some sort, the iterative updating 
of the nonlinear terms in the quasilinear form of the equations of motion 
treated by the present procedure goes a considerable way to relieving 
the restrictions imposed by the linearization inherent in other methods, 
at least in the y -direction. APPENDIX A describes the solution proce- 
dure in detail. 


The Turbulence Model 

Fully -developed turbulence model. - The fully-developed turbulence 
model originally presented by McDonald and Camarata (ref. 9) for two- 
dimensional incompressible flow, forms the basis for the transitional 
turbulence model used in the McDonald -Fish transitional boundary layer 
theory and, therefore, at this juncture it is useful to describe the 
model in some detail. The fully-developed turbulence model, which is 
described in greater detail in ref. 12, is based upon a solution of the 
turbulence kinetic energy equation. The turbulence kinetic energy 
equation is a conservation equation derived from the Navier -Stokes 
equations by writing the instantaneous quantities as a sum of mean and 
fluctuating parts. The i^ Navier -Stokes momentum conservation equa- 
tion (i =1, 2, 3, referring to the three coordinate directions) is 
multiplied by the i^* 1 component of fluctuating velocity and the average 
of the resulting three equations is taken. The three averaged equations 
are summed to obtain the turbulence kinetic energy equation. The 
derivation of the turbulence kinetic energy equation has been given by 
Favre (ref. 7) for compressible flow and approximated by Bradshaw and 
Ferris (ref. 8) to boundary layer flows ; a derivation and discussion 
of the turbulence kinetic energy equation for hypersonic flow is given 
in APPENDIX B. 



As shown in APFENDIX B, the boundary layer approximation to the 
turbulence kinetic energy equation is given by 


w (i puqZ ) + df ( T~ p v 


q 2 ) = - p uV 


advection 


"3y (PV+ -2- (pv)' 

diffusion 


l 2 ) 


du 
dy 

production 

P € 

dissipation 


- / 2 i 2 t <*u 
-/>< u 

normal stress production 


+ P 


du‘[ 


dx ,■ 

pres sure -dilitation 


(15) 


All calculations reported in the investigation were made with the usual 
assumption of zero pres sure -dilitation contribution to the energy 
balance (ref. 8) unless otherwise stated. The turbulence model is 
developed by integrating Eq. (15) with respect to y between the limits 
y = 0 and y = 8 whi ch leads to 


i lk f* P Bq2 ay= f 0 S -p u ' v ' 


du 

dy 


r 8 _ 

dy -J p e c/\j 

J ° H y (16) 




where 


z [ 2 q 2 ( P u ~ P " P 'v' + 2 ( P v) ' ^ 2 ] £ 


(17) 


Following Townsend (ref. 18) and Bradshaw and Ferris (ref. 8) structural 
coefficients a n and L are introduced, together with a mixing length/; 
these scales are defined as 


-u'v 








2 


— — 3/2 
€ - (- U V ) /L 


(-U V 


1/2 


I 


du 

dy 


(18) 
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For fully -developed turbulence the structural coefficients d| , a 2 , and 
0 3 are assumed constant having values 0.15, 0.50, and 0.20, respectively 
(refs. 9 and 19 )- As is discussed subsequently, becomes a variable 
in the transitional regime. Using Eq. (l8), Eq. ( 17 ) is put in the form 


d 

dx 


( ti Pe u e 3 

V 2 a 1 


' Pe. u e 



- 4 > 3 + 



P 


du'i 
dx j 


d y 


(19) 


where 




du/u e v 2 , 

V) 


( 20 ) 



du/ue 

drj 



( 21 ) 



I t du/u e \ 2 S du 
[IT d v ) u e dx di n 


( 22 ) 


where 77 is a nondimensional transverse distance y/ 8 *, 8 + is an arbitrary 
reference length, and 8 the boundary layer thickness. 

The left-hand side of Eq. ( 19 ) represents the streamwise rate of 
change of turbulence kinetic energy and is derived directly from the 
turbulence kinetic energy advection term. The term p e represents 
the integral of turbulence production minus dissipation and 3 

is the normal stress production. The terms designated by E are turbu- 
lent source terms resulting from disturbances imposed upon the boundary 
layer by the free stream. As shown in Eq. ( 17 ), E is the sum of two 
major contributions, the first (q 2 /2)(pud8/dx-/ov) representing the 
free-stream velocity disturbance (i.e., free - stream turbulence entrained 
by the boundary layer) and the second, pV+(/ov)' q Z /2 > representing 
the direct absorption of acoustic energy. (This acoustic energy 
absorption term will be subsequently explained in some detail.) It is 
the strength of the source term E which dictates the transition location. 
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For fully -developed turbulent flow, as in ref. 9 , L and J i are 
given by 

-g" = O.i tanh {7y/(o.l 8)] (23) 


— = tanh [Vy/loo] (24) 


where 2. CO is the "wake" value of the mixing length at any particular 
streamwise station. Although Eqs . (23) and (24) give accurate 
representations of £ and L through most of the turbulent boundary 
layer, it is well-known that they overestimate the length scales within 
the viscous sublayer and are somewhat inaccurate at low Reynolds 
numbers. Following McDonald and Fish (ref. 12) the experimentally 
observed damping effect in the viscous sublayer is modeled by assuming 
intermittent turbulence within the sublayer leading to the relation 

-77 -- r(7V) T = mao/dy)' i T ( 25 ) 


In Eq. (25) r is the intermittency factor, 2 the damping factor, and 
the subscript T indicates the value with turbulent flow. Obviously, 
is equal to the square root of T. As in ref. 12, the present 
investigation assume s th at the damping distributes normally about a 
mean height y* (y +: y-/r / p / v) with a standard deviation a leading to the 
equation 

2) = p ' /2 {(y + - y + )/cr} (26) 


where P is the normal probability function; y + is taken as 23, and cr 
as 8 . A detailed discussion of the sublayer damping treatment is 
presented in ref. 12. In the present calculations the von Karman 
constant k was taken to be 0.43. 
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In regard to the low Reynolds number effects, Coles (ref. 20) has 
observed and correlated the departure of the mean velocity profile of 
a flat plate turbulent boundary layer from the usual similarity laws 
known to hold at higher Reynolds numbers. Using Coles' correlation of 
the mean velocity profile in the low Reynolds number regime, McDonald 
(ref. 21) integrated the boundary layer equations of mean motion to 
obtain local distributions of turbulent shear stress and evaluated the 
local mixing length distributions from the assumed mean velocity dis- 
tribution and the computed shear stress distributions. Based upon these 
calculations, a low Reynolds number correction for the dissipation 
length of the form 


L = Loot I + exp (-1.63 ln R 0 + 97) ] ( 2 ?) 


was derived where Lco/S is given by Eq. (23). In the calculations 
presented in the present report the dissipation length used was obtained 
by multiplying Eq. (27) by the sublayer damping factor, 2> . 

When numerical values of the structural coefficients a n are 
specified, Eqs . (24), (26), and (27) are used to represent L and l , and 
the pressure dilitation is either neglected or modeled, the turbulence 
kinetic energy equation, Eq. (19), becomes an ordinary differential 
equation with the dependent parameter £<£>(x) which is solved in conjunc- 
tion with the boundary layer momentum and energy equations to predict 
the development of both the mean flow field and the turbulent shear 
stress. 

In addition to including the turbulence kinetic energy equation 
in the set of equations governing the boundary layer development it is 
nece ssary to specify a model for turbulent heat flux contr ibuti on, 

-~p~Cp v't 1 . As previously stated, in the present procedure, v , T / is 
specified by assuming a turbulent Pra ndtl number , Pr T , which relates 
the velocity -temperature correlation, v , T / , to the Reynolds stress, 
l/v', through Eq. (8). The turbulent Prandtl number distribution used 
in the present procedure varies with distance from the wall in the 
manner suggested by Meier and Rotta (ref. 22). At this juncture it 
should be poin ted out that an alternative procedure can be used to 
determine v't' , bas ed up on an easily derive d con servation equation for 
either the quantity T ' 2 or the correlation, v't', which is similar in 
form to the turbulence kinetic energy equation, Eq. (15). However, to 
solve this new conservation equation it is necessary to assume a uni- 
versal structure relating quantities analogous to dissipation, 
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production, etc. While sufficient experimental data exists to allow 
valid modeling of the required terms for the turbulence kinetic energy 
equation, the existing data does not indicate how proper modeling could 
be carried out for the v'T 1 conservation equation. Thus, at least for 
the present, the approach based upon a tu rbule nt Prandtl number appears 
preferable to an approach based upon the m'j' conservation equation. 

Original transitional turbulence model. - The Me Donald -Fish 
transitional turbulence model is based upon a solution of the turbulence 
kinetic energy equation with the structural coefficients modified from 
their fully-developed turbulence values. The following discussion of 
the model .condenses the presentation in ref. 12. For fully-developed 
turbulent flows the structural coefficients a, , a 2 , and d 3 are assumed 
constant and are set equal to 0.15, 0.50, and 0.20, respectively. 
Although it is probable that all coefficients vary in the transitional 
regime only the coefficient a, contributes significantly to the energy 
balance when the boundary layer is far removed from separation, as 
shown by Bradshaw (ref. 19). Therefore, the transition model assumes 
a 2 and a 3 are equal to their fully-developed values and that only q ( 
need be modified. At this juncture it should be mentioned that in a 
recent unpublished study at UARL concerned with the prediction of 
transitional separation bubbles, Briley and McDonald have found it 
necessary to modify the values for 0 2 and 0 3 slightly, however, sample 
calculations indicate that the calculations presented in the present 
report are quite insensitive to the assumed values of a 2 and a 3 . 

On the basis of current knowledge, it is expected that the 
structural coefficient a, would depend both upon the applied strain, 
dTJ/<3y and the ratio of the total effective viscosity to the kinematic 
viscosity. In regard to the first of these effects, the applied strain, 
both the experiments of Rose (ref. 23) and the analysis of Shamroth and 
Elrod (ref. 21+) indicate that for large ratios of effective to kinematic 
viscosity the adjustment of Q, from an initial value to its fully- 
strained value occurs in a few boundary layer thicknesses. Transition 
to turbulence, on the other hand, typically takes place over a distance 
of many boundary layer thicknesses and, thus, if the effects of 
viscosity are ignored, the adjustment of a ( to its fully-developed 
value under the influence of mean strain can be considered instantaneous 
in comparison with the time scale of transition. Initial attempts to 
predict transitional flows, used a fully-strained value of a t = 0.15. 
However, it soon became apparent that accurate predictions of transition 
required a much lower value for a, and, consequently, McDonald and 
Fish (ref. 12) assumed the departure of d| from its fully -developed 
value was due solely to the effect of viscosity. The relationship 
between Qj and viscosity was quantified by introducing a turbulent 
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R eynold 5 number, R r . In defining R r the velocity scale is taken as 
(-u'v 1 ) 2 and the length scale is taken as the conventional mixing length, 
L, so that 


1/2 

R r = (-U'vO Jt, 

V 


( 28 ) 


Noting that 


(-uV)= J 2 


du .2 $\j 

Jy ) = "t 77 


(29) 


the turbulence Reynolds number becomes 


Rr = 



( 30 ) 


To be consistent with the integral turbulence kinetic energy equation, 
Eq. (19)> a layer -averaged turbulence Reynolds number, R r , is introduced 
as 


R 


T 




i/dy 


( 31 ) 


where , the sublayer thickness, is defined as the location at which 
the laminar stress has fallen to four percent of the total stress (the 
four percent definition gave a sublayer mean temperature in very good 
agreement with the so-called Eckert reference temperature). 

The McDonald -Fish model assumes that the turbulence Reynolds 
number, R r , is the sole variable influencing the development of and 
a relationship between and R r is obtained by. considering the 
development of an incompressible constant pressure flat plate equilibrium 
turbulent boundary layer. It should be noted that under the assumption 
that Qj is solely dependent upon R T it is only necessary to derive a 
relation for one set of flow conditions to obtain a universally valid 
relationship. Based upon experimental observation, it is readily 
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ascertained that for the incompressible constant pressure equilibrium, 
turbulent boundary layer, , E , and d^/dx are negligible which leads 
to the reduced turbulence kinetic energy equation 


din a i din S 4a ( 4>z ® 

dR 0 ' dRg " R Q c f 4 >, 8 


(32) 


In Eq. (32), C f is the skin friction coefficient and the independent 
variable has been changed from streamwise distance, x , to momentum 
thickness Reynolds number, R^ , using the momentum integral equation. 
Equation (32) has the general form of a Bernoulli equation 


din a 
dFtg 


z f(Rff) + a, g (R Q ) 


(33) 


which has the general solution 


a,= (exp/ f dR 6 ,)/(C-/gexp/ fdR g dR g ) (34) 


When the mutual dependence between R T , f , g , and R g is determined for 
the flat plate equilibrium boundary layer, Eq. (34) provides the 
required general relationship between a ( and R r . The quantities f and 
g are evaluated by first noting that for the flat plate equilibrium 
boundary layer Q/S varies very slowly with . neglecting this varia- 
tion leads to the result 


f = 



(35) 


The quantity g which is equal to the factor -4 <p^9/( Cf<£ ( 8 ) is evaluated 
numerically by integrating the profiles of Maise and McDonald (ref. 25). 
Over a wide range of Reynolds numbers the grouping is found to be 
sensibly constant with a value approximately 6 . 66 . Thus, Eq. (34) 
becomes 
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a, = a 0 (R b /R Qq )/[i + 6.666 a 0 (R 0 /R 0q - i)] 


(36) 


where the arbitrary constant Q 0 is the value of a, when is equal to 
R$ 0< It should be pointed out that at large values of R£/R£ 0 ,a ( 
asymptotes to the fully-developed value of 0.15 • The independent 
variable of Eqs. (36) and (27) is changed from to R T by using the 
profiles of Maise and McDonald (ref. 25) which integrate to give 


R q = 68. 1 R r + 614.3 R T > 40 


(37) 


At low Reynolds numbers good results are obtained using the equation 

- 0 22 _ , 

R^ - 1 00 R f R r < l (38) 


In the intermediate range KR T <40, the two distributions, Eqs. (37) 
and ( 38 ), are joined by a cubic constructed to match the value and slope 
at the join points. Finally, the constant of integration, a Q , is 
determined on the basis of comparison between experiment and theory. 

Best agreement between theory and experiment for low Mach number, 
adiabatic wall boundary layers is obtained by setting a Q equal to 0.0115 
when R T is equal to unity. Then, from Eq. (38), =100 . 

Thus, in summary, the original transitional turbulence model is 
based upon three assumptions : 

1. The development of the transitional boundary layer Reynolds stress 
is given by the turbulence kinetic energy equation. 

2. The structural relations required to integrate the turbulence 
kinetic energy equation are identical in the fully-developed 
turbulent and transitional cases with one exception, that 
exception being the coefficient a t . 

3. The structural coefficient 0| is a universal function of the 
turbulence Reynolds number, R r . The functional relationship is 
obtained from consideration of the incompressible flat plate 
equilibrium boundary layer. 
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Modified transitional turbulence model. - McDonald and Fish 
(ref. 12) demonstrated the ability of the original transitional turbu- 
lence model to predict transitional boundary layer behavior. They 
showed that the model accurately predicts the effect of free -stream 
turbulence upon near adiabatic wall boundary layers as well as the 
length of the transition region as a function of transition Reynolds 
number in both subsonic and moderately supersonic flows. In addition, 
the theory accurately predicts mean velocity profiles, integral 
thicknesses, Stanton number, and skin friction coefficients for a wide 
variety of transitional boundary layers. Details of all these compari- 
sons are documented in ref. 12. In the present investigation the 
transitional turbulence model was used primarily in its original form, 
however, in the course of the investigation it became obvious that the 
original model did not properly predict the effect of wall temperature 
on boundary layer transition location. Although some controversy exists 
as to the effect of wall temperature on transition location, discussed 
in detail in the subsequent section, the general opinion is that 
increasing the wall temperature destabilizes the boundary layer. The 
data of Zysina-Molozhen and Kuznetsova (ref. 15 ) showing the effect of 
wall temperature upon transition location for a constant free -stream 
turbulence level in incompressible flow are presented in Fig. 1. As 
can be seen in Fig. 1, the data indicate increasing wall temperature 
destabilizes the boundary layer. In its original form the McDonald - 
Fish theory predicted heating to stabilize and cooling to destabilize 
the boundary layer, a prediction obviously in variance with the data 
and arising solely from the variation of viscosity with temperature. 

The theory was revised so that AB , the value of 2 fi/vO/q 2 at R T = 1, 
is a function of the wall -to-free -stream static temperature ratio, 

T w /T e . The functional relationship was chosen so that for ■ incompressible 
flow the theory would agree with the data of Zysina-Molozhen and 
Kuznetsova in the range 0.5< T w /T e < 2 .8. For T w /T e >2.8 the transition 
Reynolds number was assumed to be independent of wall temperature ratio. 
The variation of AB with T w /T e required to allow the theory to predict 
the foregoing effect of wall temperature upon transition Reynolds number 
for a fixed disturbance level is shown in Fig. 2. It is, of course, 
quite reasonable to hypothesize that at a given Reynolds number the 
temperature gradient effects arising from wall heating would increase 
the amount of Reynolds stress created by turbulent intensity. 
Unfortunately, the extent of such an augmentation cannot presently be 
calculated and recourse must be made to empiricism. 
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COMPARISON WITH EXPERIMENT 


Initiation of Transition 

Primarily this investigation proposes to assess the applicability 
of the McDonald -Fish transitional turbulence model in the low hyper- 
sonic Mach number regime. Necessarily, for this purpose theoretical 
predictions of transitional boundary layer behavior must be compared 
with experimental data. Although a moderate amount of data exists for 
hypersonic transitional boundary layers, much of this data is contra- 
dictory, particularly as to the effect of flow variables such as Mach 
number, wall heating, etc. upon transition location. An example of the 
difficulty encountered in obtaining hypersonic transition data is 
demonstrated by the work of Stalmach, Bertin, Pope, and McCloskey 
(ref. 26). In Run 11, Condition 4 of ref. 26, Stalmach, et al., 
investigate transition on a 12 deg cone at an edge Mach number of 7* 
Transition was determined from heat transfer measurements taken along 
two rays of the cone surface. In these apparently carefully controlled 
experiments with no evidence for significant asymmetry in the flow, 
transition location on the two rays differed by approximately 15 per- 
cent. Since the two rays would have identical conditions for the pur- 
pose of making a theoretical prediction of boundary layer development, 
the 15 percent discrepancy in the transition Reynolds number gives an 
estimate of how closely theory and data can be expected to agree. 

An important function of any transition theory is the prediction 
of transition location with various free-stream disturbances. In its 
original form, the McDonald -Fish transitional boundary layer theory 
successfully predicted the transition location as a function of free- 
stream turbulence for low-speed flows (ref. 12). The extension of this 
capability to moderately hypersonic flows was desired. In the present 
report transition is assumed to begin at that station where on a log- 
log plot straight line approximations to the Stanton number -Reynolds 
number curves in the laminar and transitional regimes intersect. 

An immediate problem arises since the theory requires a free- 
stream turbulence intensity level as input and the (NASA) data documents 
a fluctuating pressure level. For the purpose of this investigation,, 
the free-stream disturbance was assumed to be a plane wave propagating 
at a Mach number relative to the free stream. For such a disturbance 
the unsteady Bernoulli equation leads to a relation between fluctuating 
velocity and fluctuating pressure (e.g., ref. 27) of the form 
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(39) 


ui I . K | I 

"e' 7M e M r I P e I 


where, as discussed by Laufer (ref. 27), K’ 1 may be thought of as the 
integral of a space -time correlation function. The relative Mach number, 
M r , is assumed to be the function of Mqq given by Laufer (ref. 27). If 
we further assume that the disturbance is an acoustic wave, v / and u 
are related through the equation 




(UO) 


The turbulence kinetic energy is then approximated by 



v ' 2 + w ' 2 




( 41 ) 


Using the plane wave to model the disturbance, knowing the tunnel Mach 
number, Mqq, and the boundary layer edge Mach number, M e , the quoted 
free -stream fluctuating pressure, P'/P e , easily converts to an equiva- 
lent turbulence kinetic energy velocity fluctuation u'/ug which is used 
as input to the transitional boundary layer prediction procedure. Note 
that at this point it is also be ing a ssumed that the net acoustic energy 
absorbed by the boundary layer p'v' + {p\j)'q^2 in the source term 

(Eq. (17)) is negligible in comparison with free -stream turbulence 
intensity entrained by the boundary layer and thus the source term E 
of Eq. (17) is assumed to be given by 



Predictions obtained from the transitional boundary layer theory 
using both the original and modified transitional turbulence models are 
presented in Fig. 3 and compared to the data of Stainback (ref. 28). 

The original theory accurately predicts the variation of transition 
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Reynolds number with fluctuating pressure, however, this agreement 
requires a different value of K for the Mach 6 and Mach 8 tunnels. 

The Mach 6 tunnel tests and the Mach 8 tunnel test were run for 
different values of wall temperature ratio. Since the original theory 
incorrectly predicts the effect of wall cooling on transition location, 
the calculations for initiation of transition were repeated using the 
modified theory. The comparison in Fig. 3 shows that the modified 
theory compares favorably with data when the same value of K is used 
for each tunnel. The single value of K for all three tunnels obtained 
with the modified theory represents a significant improvement over the 
results obtained with the original theory since it strongly supports 
the hypothesized relation between free -stream turbulence and the measured 
fluctuating pressure. 

Rather than use free -stream turbulence intensity to trigger 
transition, it is possible to trigger transition by assuming a percent- 
age of incident acoustic energy is absorbed directly by setting q 2 at 
the outer edge of the boundary layer equal to zero and setting (pV) 
equal to a nonzero value in the source term (see Eq. (17)). For an 
incoming acoustic wave, p'v' is negative with respect to the boundary 
layer coordinate system and its magnitude can be estimated in terms of 
P'/P e from Eqs . (39) and (40). The estimate was carried out assuming 
1 and pV to be equal to |p'| jv'| (i.e., the correlation coeffi- 
cient pV/lP'l |v'| =1) • If no energy were lost as the wave passed 

through the boundary layer, the wave would be reflected with undiminished 
strength and the contrib ution of the reflect ed w ave at the edge of the 
boundary layer would be -p'v' ; thus, no net pV would result. However, 
the loss of some energy to the boundary layer could trigger transition. 
With this in mind, a study was initiated to d eter mine what percent 
acoustic energy absorption, i.e., percent of PV absorption, would 
trigger transit ion. For the purpose of this qualitative investigation 
the quantity, (p\z)'q z /2 , which is at most equal to and most likely 

less than pV > is neglected and the source term, E , of Eq. (17) becomes 

E = [ - PV j (43) 

L J e 


The results of this study are presented in Fig. 4, where the theoretical 
predictions are compared to the experimental data of Stainback (ref. 29). 
As can be seen, the theory only requires a small amount of absorption 
to trigger transition at the streamwise station where transition occurs 
experimentally (a finding in agreement with the generally accepted 
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assumption that only a small amount of acoustic energy is lost as the 
wave passes through the boundary layer and reflects from the wall). 

Thus, it appears that if the amount of acoustic energy absorbed could 
be predicted from an independent analysis, the McDonald -Fish transition 
procedure contains a mechanism for predicting transition due solely to 
absorption of acoustic energy. Needless to say, no transition theory 
by itself can indicate whether hypersonic transition is triggered by 
free -stream turbulence or absorbed acoustic energy or some combination 
of the two. However, used in conjunction with experimentally measured 
input, the theory could indicate which of the two mechanisms is most 
likely to trigger transition. 

Effect of free -stream Mach number. - The comparisons between theory 
and experiment presented in Figs. 3 and 4 were carried out at a constant 
edge Mach number, M e = 5 . An important question is, how is the 
transition location affected as the edge Mach number varies. To deter- 
mine the effect of Mach number variation on transition location it 
would be necessary to run a series of experiments in which only the 
edge Mach number varied and all other quantities, particularly the 
free -stream disturbance and the ratio of wall-to-free -stream total 
temperature, were held constant. Such experiments have been carried 
out by several investigators in the range 0<M e <4 and the data is 
presented by Zysina-Molozhen and Kuznetsova in ref. 15- According to 
the data of ref. 15, the transition Reynolds number increases with edge 
Mach number in the entire range of Mach numbers considered. A second 
set of experiments determining transition location as a function of 
fluctuating pressure level for various Mach numbers has been reported 
by Stainback, Fischer, and Wagner (ref. 29). Based upon this data and 
that of Pate and Schueler (ref. 2), Stainback, Fischer, and Wagner 
conclude that Mach number has a strong effect on transition location 
in the range 5< M e < 14 . Furthermore, the data indicates a continuously 
stabilizing effect of increasing Mach number from M e = 5 to M e = 16 at 
a constant value of fluctuating pressure ratio P'/fg . The continuously 
stabilizing effect of increasing Mach number upon transition even for 
Me<3.5 reported in ref. 15 is in contradiction with certain compila- 
tions of wind tunnel data which show that increasing the Mach number 
has a stabilizing effect for M e > 3 . 5 > but a destabilizing effect in the 
range 1.0<M e <3*5 (e.g., ref. 30). However, no reason exists to 
assume the data used in such compilations were all taken at the same 
disturbance level. In this regard, Morkovin (ref. 30) points out that 
free -flight cone data, which might be expected to have a smaller 
variation of disturbance intensity with Mach number than tunnel data, 
shows a continuously stabilizing effect of increasing Mach number even 
in the range 1.0 < M e < 3-5 - 
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A somewhat different point of view as to the effect of edge Mach 
number on transition location has been expressed by Pate and Scheuler 
(ref. 2). Pate and Scheuler assume that the dominant effect on transi- 
tion in supersonic wind tunnel experiments is the energy radiated from 
the tunnel wall boundary layers and not the Mach number per se. 

Following this assumption, Pate and Scheuler correlated transition 
location for a wide range of Mach numbers, 3 < Me < 8 5 with the wall 
boundary layer parameters and tunnel size without any explicit Mach 
number dependence. While the Pate -Scheuler correlation is quite useful 
in comparing transition experiments made in different facilities, it 
does not predict how a transitional boundary layer behaves in free 
flight, which is a problem of great interest. Furthermore, the Pate- 
Scheuler correlation is not in contradiction with a Mach number depen- 
dence of transition since disturbances emanating from the tunnel wall 
boundary layer vary with the free -stream Mach number. 

With the above consideration in mind, a series of numerical 
calculations were run using the McDonald -Fish transitional boundary 
layer theory to predict the effect of free -stream Mach number on 
transition location at a given disturbance level, uyu e = 0.03. The 
variations of skin friction coefficient through the transitional regime 
as predicted using both the original and modified transitional turbu- 
lence models for a range of Mach numbers at a free -stream turbulence 
level, u'/u e , equal to 0.03 and an adiabatic wall are presented in 
Fig. 5* Both the original and modified theories predict that increasing 
the Mach number stabilizes the boundary layer, although the modified 
theory predicts a considerably weaker stabilizing effect than that pre- 
dicted by the original theory. The predicted stabilizing effect of 
Mach number for a given free -stream disturbance level and wall tempera- 
ture ratio is believed to be in agreement with experiment, although, as 
mentioned previously, some controversy exists on interpretation of the 
various experimental data. It should be noted that in Fig. 5 the pre- 
dicted fully-turbulent value of skin friction shows agreement with the 
turbulent prediction of van Driest (ref. 31). The variation of transi- 
tion location predicted by both forms of the present prediction theory 
is compared with the data correlation over a range of free -stream 
disturbances compiled by Zysina-Molozhen and Kuznetsova (ref. 15) in 
Fig. 6. The theoretical calculations presented in Fig. 6 were carried 
out for uVu e = 0.03. The modified theory is in good agreement with 
experiment, particularly in the low hypersonic Mach number regime where 
it represents a significant improvement over the original theory. 

Effect of wall temperature . - As previously mentioned, the wall 
temperature ratio, T w /T aw , may represent another important parameter 
which influences transition location. Experimental evidence on the 
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effect of varying wall temperature is more contradictory than that of 
varying Mach number. For example, Zysina-Molozhen and Kuznetsova 
(ref. 15), van Driest and Boison (ref. 32), and Cary (ref. 33) all 
observed cooling to stabilize the boundary layer on smooth walls. On 
the other hand, Sanator, DeCar lo, and Torillo (ref. 34) observed no 
effect of wall cooling and Richards and Stollery (ref. 35 ), and Sheetz 
(ref. 36 ) observed, under certain limited conditions, a destabilizing 
effect of wall cooling. 

As in the case of the previously discussed effect of Mach number 
upon transition location, some of the contradictory evidence regarding 
the effect of wall cooling is due to variation of the free-stream 
disturbance level as the wall temperature ratio changed. However, in 
the case of wall cooling a new phenomenon may arise. As pointed out by 
van Driest and Boison (ref. 32), among others, as the wall temperature 
is lowered, the boundary layer thickness decreases making the ratio of 
the effective roughness height to the boundary layer thickness pro- 
gressively larger. As cooling continues, a wall temperature may be 
reached at which the stabilizing effect of cooling is offset by the 
destabilizing effects of roughness. In evaluating the various wall 
cooling data, it should again be emphasized that the data discussed 
by Zysina-Molozhen and Kuznetsova (ref. 15), which includes data from 
the supersonic experiments of Higgins and Pappas (ref. 37 ), specified 
the free-stream turbulence as wall temperature was varied and it was 
this data which was used to evaluate the modified transitional turbu- 
lence model. 

Several sets of calculations were performed using the McDonald -Fish 
transitional boundary layer theory to predict the effect of wall cooling 
on transition location. As shown in Fig. 7 , the original theory pre- 
dicts a destabilizing effect of wall cooling which is contradictory to 
the majority of smooth wall experimental evidence, however, the 
modified theory predicts a stabilizing effect of wall cooling. In 
regard to the predicted effect of wall cooling, it should be recalled 
that the modified transitional turbulence model was developed by 
matching theory and experiment in the incompressible regime and, thus, 
the predicted stabilizing effect of wall cooling at M g = 4.8 is not an 
automatic consequence of the model assumptions. 
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Boundary Layer Behavior Through Transition 


Prediction of Stanton number. - A major motivation in developing 
a procedure to predict the behavior of transitional boundary layers is 
the large heat flux which occurs at the end of the transitional regime. 
This large heat flux may lead to peak values of wall temperature which 
must be predicted accurately to insure structural integrity. Compari- 
sons between predictions of transitional heat transfer from the 
McDonald -Fish theory and the experimental data of Stainback (ref. 29 ) 
are shown in Fig. 8. In the case of both the 20 inch tunnel and the 
Mach 8 tunnel data the transition length scale is predicted quite 
accurately. In addition, although not shown, the theoretical predic- 
tions of transition length are in good agreement with the semiempirical 
criteria used by Harris (ref. 4) in successfully predicting a wide 
variety of transitional boundary layers. The predictions of fully- 
turbulent Stanton number for the 20 inch hypersonic tunnel presented 
in Fig. 8 are approximately ten to twenty percent below the measured 
values; however, in the laminar regime the theoretical predictions are 
in excellent agreement with the theory of Cohen (ref. 38 ), whereas the 
data is again fifteen percent higher. Since there is no apparent 
reason for the experimental heat transfer to be higher than the pre- 
dicted heat transfer upstream of transition, it was concluded that the 
20 inch tunnel experimental data was uniformly in error. The fully- 
turbulent Stanton number predictions for the Mach 8 tunnel flow condi- 
tions are in good agreement with data. It should be noted that the 
results presented in Fig. 8 differ from those of ref. 29 . The 
theoretical predictions presented in ref. 29 were preliminary predic- 
tions obtained with a theory having a slightly different sublayer 
model than that used in the present calculations. A second set of 
Stanton number predictions carried out for the data of Holloway and 
Sterrett (ref. 39) are presented in Fig. 9* For this comparison pre- 
dictions of both the length of the transition region and the Stanton 
number distribution itself are in good agreement with data. In con- 
trast to the case of the Stainback data, although a discrepancy between 
theory and experiment does appear in the peak value of Stanton number, 
the predicted maximum Stanton number in this instance exceeds the 
experimental value by approximately 15 percent. While no definite 
conclusion can be reached concerning the accuracy of the results in 
the fully-turbulent regime, the theoretical predictions are in good 
agreement with the flat plate analysis of van Driest (ref. 31); whereas 
the experimentally observed values are lower. For the above predic- 
tions of transitional Stanton numbers the input free -stream turbulence 
level was chosen to give agreement between the theoretically predicted 
and experimentally observed transition locations. This was also done 
in obtaining the predicted transitional velocity profiles subsequently 
presented. 
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Prediction of mean velocity profiles. - A stringent test for 
assessing the ability of a boundary layer calculation procedure to 
predict transitional behavior is the comparison between predicted and 
experimental mean velocity profiles. For this reason three comparisons 
between measured and predicted velocity profiles were made during this 
investigation. The first of these compares the theoretical predictions 
of the McDonald -Fish transitional boundary layer theory with the 
experimental data of O'Donnell (ref. 40). For this comparison, as for 
all other comparisons where no free -stream turbulence was. quoted, a 
turbulence level of u'/u e was chosen which initiated transition at the 
experimentally observed streamwise location. 

A comparison between the theoretically predicted and experimentally 
measured velocity profiles is presented in Fig. 10. Some minor 
discrepancies between theory and experiment exist in the profile wall 
region where measurement difficulties would be expected, however, 
elsewhere the predictions are in excellent agreement with the laminar 
theory presented by O'Donnell (ref. 40) and the discrepancy in the 
transitional regime does not exceed that in the laminar regime. 
Therefore, the theoretical prediction should be regarded as good. A 
comparison between theoretically predicted and experimentally measured 
momentum thickness is presented in Fig. 11 which again shows good 
agreement between theory and experiment. 

A second set of calculations was carried out to compare the 
McDonald -Fish theory with the Michel and Schmitt profile data (ref. 4l) 
presented in Fig. 12. Both the fully -laminar and fully -turbulent 
predicted profiles are in good agreement with data, as is the length of 
transition region, however, within the transitional regime itself the 
agreement is not very good. The data was measured on a cylindrical 
rod with a radius only four times the boundary layer thickness. No 
account of the finite value ratio of the rod radius to boundary layer 
thickness was taken in the calculation and this may explain some of 
the discrepancy. A slight adverse pressure gradient in the experiment 
represents a second possible source of error. The experimentally 
determined edge Mach number distribution showed considerable scatter 
between M e = 6.65 and M e = 6.55 and for the purpose of the calculation 
the edge Mach number was assumed constant at a value of 6.6. Compari- 
son of integral thicknesses and Stanton number are presented in Figs. 

13 and l4. As can be seen in Fig. l4, the measured Stanton numbers 
do not agree with either the McDonald -Fish theory or the theory of 
Cohen in the laminar regime, thus making the transitional comparisons 
doubtful. 
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The final comparison between theoretically predicted and 
experimentally measured profiles was made for the NASA Langley Research 
Center data of Fischer and Maddalon (ref. 42), which contains both 
mean velocity and density profiles. Comparisons between the theory and 
experiment for the velocity profiles is presented in Fig. 15 and for 
the density profiles in Fig. l6. Both the predicted length of the 
transition region and the velocity profiles are in reasonably good 
agreement with experiment. Although the density profiles are not as 
well predicted as the mean velocity profiles, they are still quite 
acceptable, particularly considering the large density gradients 
through the midportion of the boundary layer. General agreement is 
shown between measured and predicted displacement thickness presented 
in Fig. 17* In summary, it seems that considering the uncertainties 
in the measured data, the comparisons between the predicted and measured 
mean velocity profiles are quite good. 


The Role of Fluctuating Pressure 


Two major assumptions which are made in deriving the turbulence 
kinetic energy equation should be reexamined in the hypersonic flow 
regime. These concern the direct roles of fluctuating density and 
fluctuating pressure. As shown in APPENDIX B, the most likely direct 
contribution of fluctuating density to the turbulence kinetic energy 
b alan ce occurs i n the turbulence production mechanism where the terms 
vp 'u' dU/dy and p'u'v 1 du/dy are neglected in comparison to pu 1 v'd u/d y. 
The ratios of the neglected to the included terms are given by 


_V_ P^_ 

P u V 



p' uV 
p u' v' 



(44) 


In the boundary layers investigated in the present investigation it - 
appears likely that these ratios remain small and, thus, it' does not 
appear that in the low hypersonic Mach number regime fluctuating' 
density makes a major contribution to the turbulence energy balance. 

The role of fluctuating pressure terms, however, is not as 

obvious. As previously discussed, pressure fluctuations may play a 

major role in the transition process by acting as a source term in the 

turbulence kinetic energy equation.' The calculations presented in 

Fig. 4, show that only approximately 0.3 percent of the incident 

acoustic energy need be absorbed by the boundary layer in order to 

trigger transition at the observed location in absence of any free- 

2 

stream turbulence kinetic energy, q . A second mechanism through which 
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pressure fluctuations ma y influence transition is represented by the 
pressure dilitation term P x d u' j /dx -. Although this term is usually 
neglected in forming the turbulence energy balance, there is no reason 
to expect it to be insignificant in high Mach number flows. An 
estimate of its magnitude is obtained from the continuity equation in 
the following manner . When divided into mean and fluctuating parts 
the mean flow continuity equation is written as 

IT ( P + P‘) + { P + P " 1 (Q i + u i )r + 

dpu- t dp'uj aouj dp'u\ 

4" , ■+■ ~ -f* 7 + - - 0 

dx\ dx,- dx ,• dxj 


dp 1 

at 


(45) 


Averaging Eq. (45) yields 


dp- <^Ui+/>Vi 

at + Sxj 


(46) 


which when subtracted from Eq. (45) leads to the result 
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As pointed out by Bradshaw and Ferris (ref. 8), Taylor's hypothesis 
leads to the conclusion 
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and, therefore, 


p/ dui . ^pL jtij_ _ P^J_ JLL + P> d P ' u ' i (50) 

<3x j p (3xj p dx i dx\ 


In view of the recent data of Owen and Horstman (ref. 43), Taylor's 
hypothesis must be used with caution in the transitional regime, 
however, the hypothesis is still good enough for the present order of 
magnitude arguments to be valid. When third order correlations and 
the mean flow dilitation is assumed small and the boundary layer 
approximations are made, the final expression for the pressure dilita- 
tion term becomes 


/ d u i _ _ PV d£_ 

dx j p dy 


(51) 


If a suitable model for P / v / is hypothesized, the approximation to the 
pressure -dilitation term can be included in the turbulence kinetic 
energy equation. 


In hypothesizing a model for the factor P / v / to be used in the 
present investigation to determine the contribution of the pressure - 
dilitation term in the turb ulence kinetic energy equation, the assump- 
tion was made that the P' v 7 to be used should be that generated in the 
boundary layer rather than that which propagates in from the free 
stream and is absorbed by the boundary layer. Based upon the data 
obtained by Bradshaw (ref. 19) and Laufer (ref. 44), p' v ' was as sume d 
constant across the boundary layer. The assumption of constant P 7 v' 
with distance from the wall is obviously in error a t th e wall where 
P ’ v 7 must be zero, however, the data doe s indicate P 7 v 7 is constant in 
the outer region. The magnitude of p'v was obtained by noting that 
Kistler and Chen (ref. 45) measured | P w j /r w in fully -turbulent flow to 
be a function of edge Mach number, M g . For transitional flow this 
relation was modified to 


P 7 = 


w 


(M e ) 


,8 . 

/ o S V T + vdsf 


(52) 
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where v is the kinematic viscosity, v T is the turbulent kinematic 
viscosity and f(M e ) is obtained from the Kistler-Chen data. 

The fluctuating velocity v' was estimated from 


(53) 

2 

which is known to overestimate the value of v' in trans iti onal flows. 
Using these assumptions and assuming, furthermore, that pV / I P '\ I v ( l = 
1.0 the normalized pressure -dilitation term becomes 

eL iHi . IL f5 4) 

Pe u i <5xi ‘ > e u e 3 P d y w ’ 

where | P | and (v'j are obtained from Eqs. (52) and (53), respectively. 
Since P'v' is positive (ref. 19), the pressure-dilitation term acts as 
a turbulent sink. A comparison of calculations carried, out, including 
and omitting the pressure-dilitation term, is presented in Fig. 18. 

The calculations presented in Fig. 18 were made assuming the turbulent 
source term, E, of Eq. (17) to be given by 



(55) 


2 

The magnitude of q was chosen by relating the fluctuating velocity to 
the measured fluctuating pressure through Eq. (39) • It should be noted 
that Eq. (39) contains an empirical constant which was evaluated by 
comparing theoretical -calculations which neglect pressure-dilitation 
effects with experimental data. Since the data presented in Fig. l8 
was one of the sets of data used in evaluating the empirical constant, 
the theory neglecting pressure-dilitation effects shows better agree- 
ment with data than the theory including pressure-dilitation effects. 
However, the important conclusion to be drawn from Fig. l8 is that 
although.it was given its maximum reasonable value, inclusion of the 
pressure-dilitation term made only a minor difference in transition 
location and no difference in the shape of the Stanton number -Reynolds 
number curve. Thus, at the present time in view of the uncertainty 
of the assumptions which must be made in modeling the pressure-dilitation 
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effect, little advantage accrues to the inclusion of the pressure- 
dilitation term in the turbulent energy balance for these conditions 
where M g = 5.0 and T w /T 0 = 0.6. 


The Precursor Effect 

• Hot wire anemometer studies show that significant oscillations 
occur in the outer portion of the boundary layer far upstream of the 
transition location determined by wall instrumentation (refs. 29, 46 
and 47). These oscillations, termed the 'precursor effect', are 
associated with the onset of transition, and initially appear in the 
vicinity of the critical layer and spread at an angle of approximately 
0.5 deg relative to the wall. Thus it is expected that the precursor 
effect would first be noticed in hypersonic flow at approximately 150 
boundary layer displacement thicknesses upstream of the streamwise 
station at which transition is recognized from heat transfer measure- 
ments. The McDonald -Fish theory in solving the turbulence kinetic 
energy equation in nominally laminar flow also predicts such a precursor 
effect. According to the theory a small amount of Reynolds stress 
first appears in the outer portion of the boundary layer far upstream 
of the station at which the wall properties depart from their laminar 
variation. The region of significant Reynolds stress spreads as the 
flow proceeds downstream and when this region approaches the wall, the 
wall properties such as Stanton number and skin friction differ from 
laminar. A plot of computed normalized stress profiles for conditions 
corresponding to those for Stainback's Mach 8 tunnel data (ref. 29) at 
several selected stations, presented in Fig. 19, shows that at 200 
displacement thicknesses upstream of the transition location a small 
amount of kinematic Reynolds stress appear? in the outer portion of the 
boundary layer. The transition location, x T , was again determined by 
the intersection of straight line fairing in a log-log plot of Stanton 
number versus Reynolds number. As the flow proceeds downstream, the 
Reynolds stress increases in mangitude and transverse extent. Signifi- 
cantly the predicted transverse location at which the precursor effect 
first appears is in the outer portion of the layer in qualitative agree- 
ment with experimental data. Also, the predicted spreading angle for 
the particular case considered is approximately 0.3 deg, also in 
qualitative agreement with experimental data (ref. 47). 
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DISCUSSION 


The problem of predicting the initiation of transition and the 
behavior of boundary layers through the transitional regime is one of 
the most difficult problems encountered in boundary layer theory. 
However, the contradictory experimental evidence as to the effect of 
Mach number (e.g., refs. 2, 15, and 30) or wall cooling (e.g., refs. 15? 
30, and 32 through 36) clearly indicate the need for developing a tran- 
sition prediction procedure based upon as low a degree of empiricism 
as possible. In addition to the effect of Mach number and wall cooling 
other phenomenon not considered in the present investigation such as 
strearawise pressure gradient, roughness (e.g., refs. 30, 39> and 48 
through 50) and unit Reynolds number (e.g., refs. 2, 3, 29, 30, and 48) 
may also have important effect on transitional boundary layer behavior. 
Tunnel effects which further complicate the data assessment problem 
(refs. 2, 3j and 29 ) may make valid comparisons between data taken in 
different wind tunnels difficult and, more importantly, may make the 
correct application of wind tunnel data to free -flight conditions doubt- 
ful unless the wind tunnel data is obtained under stringent control. 

With these considerations in mind, the advantages of a general 
transitional theory such as that described here are evident. In the 
present report, the theory has been examined in four respects: the 

ability to predict (1) the initiation of transition; (2) the effect of 
parametric variation upon transition; (3) heat transfer through tran- 
sition; and (4) velocity and density profiles through transition.' The 
calculations to ascertain if the procedure could predict transition 
location as a function of the free -stream disturbance level were based 
upon the free-stream disturbance as a vorticity fluctuation, i.e., 
free-stream turbulence intensity ^ At higher Mach numbers , M e > 2 .5 . , 
in addition to the free-stream turbulence the disturbance field con- 
tains both a sound mode (refs. 2 and 3) and an entropy mode. Although 
the entropy mode is unlikely to be important except at very high Mach 
numbers, it is possible that the sound mode may be the principal tran- 
sition trigger, particularly with small vorticity fluctuations. How- 
ever, even a dominant acoustic mode does not preclude free-stream 
fluctuating velocity entrainment by the boundary layer as the source of 
transition, sinca the fluctuating velocity field associated with the 
acoustic field could be entrained by the boundary layer and cause tran- 
sition in the same manner as free-stream turbulence intensity causes 
transition at low Mach numbers. The calculations presented in Fig. 3 5 
which are based on this concept of an associated ’ fluctuating velocity 
field, are encouraging since the same value of K gives good agreement 
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between theory and experiment for three tunnels. Obviously, the theory 
would have to be tested over a much wider range of experiments before 
any definite conclusions could be made. 

An alternative and direct source for transition is the absorption 
of the incident acoustic energy by the boundary layer. At present no 
way exists to estimate accurately what fraction of the energy of the 
incident acoustic wave is absorbed by the boundary layer, although this 
fraction absorbed must be small. The McDonald -Fish transitional theory 
predicts that less t han one percent of the acoustic energy; i.e., less 
than one percent of P ' v', need be absorbed to trigger transition (see 
Fig. 4). However, until the development of an adequate theory to pre- 
dict acoustic energy absorption by a boundary layer, it is impossible 
to judge the reasonableness of a one percent energy absorption. In 
regard to the ability of the McDonald -Fish theory to predict the effect 
of parametric variation upon transition, predictions made with the 
modified theory for the variation of transition location with Mach 
number agree with the data compiled by Zysina-Molozhen and Kuznetsova 
(ref. 15), as shown in Fig. 6. In addition, the modified theory 
correctly predicts the effects of wall cooling in the absence of rough- 
ness, Fig. 7, at high Mach numbers. In regard to the other parameters, 
calculations presented by McDonald and Fish (ref. 12) show the theory 
successfully applied in the presence of streamwise pressure gradients 
and wall roughness for low speed flow. Based upon the success of the 
present investigation, these trends could probably be predicted 
accurately in the low hypersonic Mach number regime . The ability to 
predict variation of transition location with changes in flow para- 
meters from the fundamental conservation equations is one of the chief 
advantages of the McDonald -Fish transitional theory. 

Predictions of the Stanton number distribution through transitional 
flow indicate the theory accurately predicts the length of the transi- 
tional region in the low hypersonic Mach number regimes. Predictions 
of the Stanton number distribution itself appear to be good. In the 
case of the Stainback data (ref. 29), as shown in Fig. 8, within 
limits of experimental scatter, the theory predicts both the length of 
transition and the Stanton number accurately. 

The final test of the McDonald -Fish transitional theory is the 
ability to predict mean velocity profiles. Although the predictions 
shown in Fig. 12 were not in good agreement with the data of Michel 
and Schmitt (ref. 4l), the predictions shown in Figs. 10, 15, and 16 
agreed with the data of O'Donnell (ref. 40) and the NASA Langley data 
of Fischer and Maddalon (ref. 42). The velocity profiles were calcu- 
lated using the original theory, but no significant changes in this 
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regard are expected from the modified theory. Thus, on balance it is 
concluded that t'he McDonald -Fish transitional -theory can accurately 
predict: (l) the length of transition, (2) transitional heat transfer, 

(3) transitional velocity profiles, and (4) effect of Mach number, wall 
temperature, and free -stream disturbance level upon relative transition 
location. In addition, the results are encouraging as to the theory’s 
ability to predict quantitatively, the absolute transition Reynolds 
number corresponding to a specified set of conditions and the theory 
predicts a precursor effect in qualitative agreement with that experi- 
mentally observed. 


CONCLUSIONS 


1. The McDonald -Fish transitional boundary layer theory shows 
considerable promise for predicting the transition Reynolds number at a 
given free -stream disturbance level in the low hypersonic Mach number 
regime . 


2. The theory accurately predicts the change in transition loca- 
tion caused by the free -stream Mach number, wall temperature ratio, and 
free -stream disturbance level. 

3. The theory contains mechanisms for triggering transition by 
free -stream turbulence entrainment or from direct acoustic energy 
absorption. 

4. The theory accurately predicts the variation of mean velocity 
profiles, Stanton numbers, and integral thicknesses through the transi- 
tional regime for low hypersonic Mach number boundary layers. 

5- The direct effects of fluctuating pressure through the 
pressure -dilitation contribution to the turbulence energy balance do 
not significantly alter the transition Reynolds number for low hyper- 
sonic Mach number boundary layers. 

6. The theory predicts turbulence production in the boundary 
layer far upstream of the place where it is first noticeable at the 
wall (the precursor effect) similar to effects noticed experimentally. 
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RECOMMENDATIONS FOR FUTURE WORK 


Further development and assessment of theoretical approaches for 
predicting transitional boundary layers must rely heavily on experi- 
mental guidance. New information of the type not previously generally 
measured in the transitional boundary layer, such as the development of 
the Reynolds stresses and the precise makeup of disturbance modes is 
urgently required. Given the very close relationship observed at low 
speeds between the intensity of the disturbance and the transition 
location, it is very surprising to find transition studies at high 
speeds where not even the free -stream disturbance level is measured. 

Insofar as the present transitional theory is concerned, should 
further improvement in the predictions be found necessary, three 
possible areas of improvement can be investigated. The first of these 
concerns the turbulence structural hypothesis assumed for the transi- 
tional boundary layer and here more refined turbulence relationships 
can obviously be examined. Secondly, if sufficient confidence can be 
established in a pressure -velocity diffusion relationship, then it may 
prove rewarding to replace the integral representation of the turbulence 
kinetic energy equation by a finite -difference solution to the actual 
local partial differential equation. However, it is felt that the use 
of an integral form of the turbulence kinetic energy equation at the 
present time avoids the difficulty of explicitly specifying this 
largely unknown pressure -velocity relationship at the same time as 
simplifying the computation. These considerations lead to what is felt 
to be a more rewarding third area of investigation; that is, of replacing 
the turbulence kinetic energy equation used in the present theory by 
the three component intensity equations and the Reynolds apparent shear 
stress equation, again all initially in integral form. A preliminary 
investigation of this four equation integral representation has shown 
encouraging results in fully -turbulent boundary layer flow. 
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APPENDIX A 


The Calculation Procedure 

The procedure used to solve the momentum and energy equations is 
a Hartree-Womersley calculation procedure similar to that described by 
Smith and Clutter (ref. 51 )• In the calculation procedure the 
differential equations are transformed to a form more convenient for 
computer solution by introducing new variables 


V r y/8 + 

(A-l) 

P u 

(A -2) 

Pe u e 


T e °-f° 


^e 0 - t ref 

(A -3) 

p ■■ p„/p 

(A -4) 


where the primes indicate differentiation with respect to 77 , T °is the 
stagnation temperature, T REF is a specified constant reference tempera- 
ture, and 8 + is a length scale usually taken to be the displacement 
thickness. In certain flows, such as the boundary layer developing on 
the cold wall nozzle, the displacement thickness may become small or 
even negative and in these cases 8 + is taken to be a linear combination 
of the displacement thickness and a constant reference length such 
that 8* remains positive. 

The momentum equation, Eq. ( 9 ), and the energy equation, Eq. (10), 
are solved by first eliminating pM through the continuity equation, 

Eq. (3), and replacing the variables T °, p, and u by G 1 , F, and /3. The 
streamwise derivatives are then replaced by finite differences leading 
to equations of the form 

A 3 F'"+ A 2 F"+ a, F'+A 0 F:A 4 (A-5) 
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B 3 O'"* B 2 g" + B| G'= B 4 


(a-6) 


where A n and B n are functions of F, G ; and their derivatives. The 
equations are linearized by assuming values for F, G 7 and their deriva- 
tives based upon the solution at the previous streamwise stations and 
the resulting linear differential equations are solved by Gaussian 
.elimination using the boundary conditions 


F',0) = , . (A-7) 

re u e 

F'<8)=0 (A -8) 

r 0 r 

G'(0) = _ W OR- G"(0) = 0 (a -9) 

T e _ 'ref 

G'(S th ) : 0 (A-10) 


where 8 TH is the thickness of the thermal boundary layer. In the 
calculations presented in the present report, u w was always taken as 
zero. Having obtained the solution of the linearized equations, the 
distributions of F and G 1 are compared to the distributions used to 
evaluate the nonlinear coefficients, A p and B p . If the old and new 
distributions agree to within a specified tolerance, the procedure 
moves to the next streamwise station. If the two do not agree, the 
procedure is repeated. In the case of turbulent flow, the coefficients 
A n and B n , also depend on the turbulent kinematic viscosity, Vy , and a 
turbulence kinetic energy equation is included in the iteration loop. 
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APPENDIX B 


The Turbulence Kinetic Energy Equation 

The turbulence kinetic energy equation is derived from the 
Navier-Stokes equations by multiplying the i"*"* 1 Navier-Stokes momentum 
equation by the i^B component of fluctuating velocity and averaging 
the results (ref. 8). If the three average equations are then summed 
the result is 



+ (pufuj + uj />X 


+ p'u'u'jj 


du\ 
dx j 

(B-l) 


+ + 1/2 



+ 



du'i 
dx j 


/_3ui 

8x| 


= 0 


where q 2 is the turbulence intensity defined by 

2 ,2 , 2,2 , . 

q-u +v +w (B-2) 


and 


ct - is the fluctuating viscous stress tensor given by 


■j (3x k «i n 


du\ 


dx y dXj 


(B-3) 


where X is the second coefficient of viscosity and 8.. is the 
Kronecker delta defined as 


Sjj=0 i ^ j 

S ii = I i = ] 


(B-4) 
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Since it is assumed that 


« I 


(B-5) 


u 

rT « I 
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the first term is approximated as 


H + p‘ a i) 




dx . 


/DU 


dx 


+ + P'*' ) 


Similarly the next set of terms is expanded as 
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Since by the boundary layer approximations 
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Equation (B-8) becomes 
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From the boundary layer approximations 
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and from experimental evidence 


p o' v' « p u'o' » p vV 


(B-13) 


Although only a small number of measurements for p'/p exist in the low 
hypersonic Mach number regime, Harvey, Bushnell, and Beckwith (Ref. 52) 
have shown that in this regime the available evidence indicates 


o < 0.2 

P 


(B-14) 


Eq. (B-l4) is assumed valid for flow conditions considered in the 
present report. Thus, 


(p u / u j + u j p' u i + p' u - /u * P uV r 


, — ~r~i du . — ~ 7 d v , - 7 — <3u - — r^r 3u 
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Although the second and third terms are usually small compared to the 
first, they are retained since they represent production by normal 
stresses in incompressible flow and may become important in the region 
of separation or transition. For subsonic and moderately supersonic 
flow in the absence of wall heating, order of magnitude arguments 
indicate the last two terms are small compared to the first term 
(ref. 8) and they are neglected in the present formulation, (in the 
case of high Mach number flows or flows with significant wall heat 
transfer, it may become important, if I p I/p ~| v '\/ v. ) With the above 
approximations, Eq. (B-ll) becomes 


\ <3U; 

V°/ + v' u l + /AS/^rr 


+ p uV 


, - I / 

+ p V V 


Following ref. 8 it is assumed that viscous diffusion is negligible. 
Thus the turbulence kinetic energy equation becomes 



+ l/2(/>v)V 


It a- 1 


Using the equation of continuity the turbulence kinetic energy equation 
is put in final form as 



advection (mean flow convection) production pressure-dilitation 

(B-18) 

- ^r( p ' v ' + 1/2 (p^'q 2 ) -pe -p(u'-v' 2 ) 

diffusion dissipation normal stress term 
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where e represents the sum of the turbulent dissipation terms. In 
obtaining Eq. (B-l8) the normal stress term (which did not make a 
noticeable contribution to the turbulence kinetic energy balance for 
calculations presented in the present report) it was assumed that 


8u 

dx 



(B-19) 
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Figure 2- - Variation of structural coefficient with temperature ratio. 
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Figure 3. - Effect of free-stream fluctuating pressure on boundary layer 
transition location on sharp cones at M e = 5. 
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Figure 4- - Comparison between measured and predicted transitional 
heat transfer with transition triggered from pressure - 
velocity fluctuation. 
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Figure 5. - Predicted variation of transitional skin friction for various Mach numbers. 






Figure 6. - Variation of transition Reynolds number with Mach number 
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Figure 7. - Predicted variation of transitional skin friction for various wall temperature ratios. 



STANTON NUMBER, St 


. LAMINAR THEORY (REF. 38) 

— MODIFIED UARL THEORY 

ORIGINAL UARL THEORY 

STAINBACK'S DATA (REF. 29) 

20 IN. HYPERSONIC TUNNEL MACH 8 TUNNEL 

M M =6.0 Moo = 8.0 

° st Laminar ^ s, laminar 



10 6 2 5 10 ? 2 

REYNOLDS NUMBER, Re 

' X 


Figure 8. - Comparison between measured and predicted transitional 
heat transfer at Me = 5. 
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Figure 10. - Comparison between measured and predicted transitional velocity profiles at Me = 24 
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Figure 11. 


Comparison between measured and predicted 
momentum thickness at M e - 2.4. 
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ORIGINAL UARL THEORY 



Figure 12. - Comparison between measured and predicted transitional velocity profiles at M 
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Figure 13. - Comparison between measured and predicted variation of 
boundary layer integral parameters through transition at 
M e - 6.6. 
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Figure 


14. - Comparison 
transitional 


between measured and predicted 
heat transfer at M e -6.6. 
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Figure 15. - Comparison between measured and predicted transitional velocity profiles at M e -6.2. 
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Figure 19. - Reynolds stress profiles upstream of transition. 
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